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We present a class of systems for which the signal-to-noise ratio as a function of the noise level 
may display a multiplicity of maxima. This phenomenon, referred to as stochastic multiresonance, 
indicates the possibility that periodic signals may be enhanced at multiple values of the noise level, 
instead of at a single value which has occurred in systems considered up to now in the framework 
of stochastic resonance. 
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PACS numbers: 05.40. +j 

In recent years, the phenomenon of stochastic reso- 
nance (SR) has been subject of intense activity, to the 
extent that many examples have been found in different 
scientific areas |i|-|l5|]. The phenomenon has been char- 
acterized by the appearance of a maximum in the output 
signal-to-noise ratio (SNR) at a non-zero noise level. In 
this sense, noise plays a constructive role since an opti- 
mized amount is responsible for the enhancement of the 
response of the system to a periodic signal, which other- 
wise would be manifested with more difficulty In spite 
of the efforts devoted to its understanding, there is an 
aspect which has not been considered up to now. Under 
some circumstances, the SNR may display a multiplicity 
of maxima and hence, there is a set of values of the noise 
level at which the response of the system is enhanced. 
In this Letter we address precisely this possibility for the 
appearance of those maxima, then we show that the con- 
cept of SR is more general than the one we already know. 
In this regard, we have found a class of systems whose 
output SNR is a non-trivial periodic function of the log- 
arithm of the noise level. Similarly, our analysis also 
reveals that the SNR as a function of the noise level may 
exhibit any number of maxima. 

Consider systems with only one relevant degree of 
freedom whose dynamics is described by the following 
Langevin equation 



dx 



-f(x,t)x + VD£(t) 



(1) 



where f(x,t) is a given function, £(t) is Gaussian white 
noise with zero mean and second moment = 
5(t—t'), and D is a constant defining the noise level. Here 
the input signal enters the system through f(x,t), and 
we will assume it to be periodic in time with frequency 
luq. The output of the system is given by v(x) = \x\ n , 
with n a positive constant. The SNR is defined, as usual, 
by 



SNR = S{ujo) v ( x) /N(ujo) v{x) 



(2) 



where S(ujq) v ^ and N(u>o) v ^ are the output signal and 
output noise corresponding to v(x), respectively. 

Under the transformation x — e^x and D — e 2l D, 
with 7 constant, Eqs. ([!]) and (Q) remain unchanged if 



f(x,t) = f(xe^,t) 



(3) 



Consequently, for the class of systems for which the pre- 
vious equality holds, for a certain value of 7, the SNR 
has the same value at D and at e 2l D. This fact occurs 
when f(x,t) = q(ln(x),t), where q is a periodic function 
of its first argument, with periodicity 7 if 7 is the lower 
positive number satisfying Eq. (||). Therefore, the SNR 
is a periodic function of the logarithm of the noise level. 
Notice that both the signal and noise are not invariant 
under this transformation, but they are changed in the 
following fashion 



— „27n 



(4) 



In order to illustrate the previous results we have ana- 
lyzed some representative explicit expressions of f(x,t). 
To this purpose we have numerically integrated the cor- 
responding Langevin equation by means of a standard 
second-order Runge-Kutta method for stochastic differ- 
ential equations Eg]. Moreover, as the output of the 



system we have used v(x) 
the case in which 



We will first consider 



f(x,t) = @ T (log 10 (x 2 ))[f3 + acos(u} t)] 



(5) 



where a and ujq are constants, and 0t(s) is a square 
wave of period T defined by 



e T (s) 



fci if sin(27rs/T) > , 
fc 2 if sin(27rs/T) < , 



(6) 



with k\ and ki constants. In Fig. 0(a) we have plotted 
the SNR corresponding to the previous form of f{x, t) , for 
particular values of the parameters. This figure clearly 
manifests the periodicity of the SNR as a function of the 
noise level and the presence of multiple maxima at D = 
Doe mT , with m being any integer number and Do the 
noise level corresponding to the maximum with m = 0. 
We show both the signal and noise in Fig. [l](b). This 
figure also corroborates the dependence of the signal and 
noise on D given in Eq. (^). 

From this example one can infer the mechanism re- 
sponsible for the appearance of this phenomenon. Due 
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to the fact that the SNR has dimensions of the inverse 
of time |L5] , its behavior is closely related to the charac- 
teristic temporal scales of the system. Thus variations of 
the relaxation time manifest in the SNR. In this example, 
when T is sufficiently large, for some values of the noise 
level the system may be approximated by 

dx 

— = -ki[/3 + asm(u t)}x + V2D£(t) , (7) 

where i = 1,2, depending on the noise level. In such a 
situation the SNR is given by 

SNR = g(a,u K 1 )k i , (8) 

with g a dimensionless function [ |l5|] . For a sufficiently 
low frequency, the SNR is proportional to ki (SNR rj 
g(a,0)ki), i.e. proportional to the inverse of the re- 
laxation time. Consequently, there are two set of val- 
ues of D for which the SNR differs in approximately 
101og 10 (fci/&:2) dB, as one can see in Fig. 0(c). We then 
conclude that multiple maxima in the SNR appears as a 
consequence of the form in which the relaxation time of 
the system changes with the noise level. 

Let us now consider another explicit expression for 
f(x,t) that mainly differs from the previous one in the 
form in which the input signal enters the equation, 

f(x,t) = fcsin(27rlog 10 (x 2 )/T) + [/3 + a cos(ui t)] . (9) 

Here the spatial and temporal dependence of f{x,t) ap- 
pears in an additive fashion instead of in a multiplicative 
way as in Eq. (^J). The results for the SNR are shown in 
Fig. (||) and also corroborate, in this case, the periodic 
dependence on the logarithm of the noise level. 

The importance of the class of systems discussed pre- 
viously lies in the fact that the periodicity of the SNR 
can be shown analytically by simple considerations about 
the invariance of the system under stretching transfor- 
mations. Numerical analyses have also revealed that the 
presence of multiple maxima in the SNR as a function of 
the noise level is a more general phenomenon than the 
situation described by the equality Thus, the ap- 
pearance of multiple maxima is robust upon variations 
of the form of f(x,t). For instance, if the periodicity in 
s of the function q(s, t) is lost for high and low values of 
s then the SNR may be a periodic function of ln(_D) for 
a bounded range of values of D. In this regard, we have 
performed numerical simulations and we have observed 
that it is possible to obtain any number of maxima de- 
pending on the form of f(x,t), even if the periodicity of 
q(s,t) does not hold for any interval of s. As an illustra- 
tive example we will analyze a case in which the input 
signal enters f(x,t) in a additive fashion as well as in Eq. 

®. 

We will consider that the form of f(x, t) consists of two 
contributions; namely, one that comes from a time depen- 
dent parabolic potential and the other without temporal 
dependence. Then, 



f(x,t) = f (x) + [(3 + a cos(cjt)} . (10) 

Parabolic potentials may arise in many physical situa- 
tions of interest. For instance, around a minimum most 
of the potentials may be approximated by a parabolic 
one. Experiments of a temporal variation of the inten- 
sity of the potential are common, this is the case of a 
single dipolc under the influence of an external oscillat- 
ing field p"7j ] . Another remarkable situation described in 
the same way corresponds to some systems around a bi- 
furcation whose control parameter varies periodically in 
time, as for example Rayleigh-Benard convection when 
the temperature difference between plates varies slowly 
in a periodic fashion |l^,[l{|. The term fo(x) may then 
represent a perturbation to this ideal situation described 
by only a time-dependent harmonic force. 

To be explicit and for the sake of simplicity, we will 
consider fo(x) to be a linear piecewise function defined 
by 
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In Fig. |3|(a) we have depicted the potential Vq(x) corre- 
sponding to the force —fo(x)x for particular values of the 
parameters ki and a. It is worth emphasizing that we 
have considered other forms of fo(x) but the same qual- 
itative results are obtained provided that fo(x) exhibits 
the main characteristics of Eq. (|Tl|), i.e. the potentials 
must have a similar location of maxima and minima. In 
Fig. ||(b) we have represented the SNR corresponding to 
the previous function fo(x) as a function of the logarithm 
of the noise level. This figure clearly displays two max- 
ima, one more pronounced than the other. Depending on 
the values of the parameters the second maximum may 
disappear [Fig. ^(c)] or it can become more pronounced 
[Fig. 1(d)]. 

In summary, we have found, for the first time, a class 
of systems for which the response to a periodic force is 
enhanced, not only with the addition of an optimized 
amount of noise, but also at multiple values of the noise 
level. Thus, the SNR exhibits a series of maxima dis- 
tributed periodically, as a function of the logarithm of the 
noise level. This feature has been found to be even more 
general since any number of maxima may be present. 
Among others, an applied aspect to be emphasized con- 
cerns the possibility for the design of devices for which 
the enhancement of an external signal may occur at dif- 
ferent values of the noise and not only at one particular 
value. Our findings, then, contribute to a wider under- 
standing of the phenomenon of SR by extending its scope 
and perspectives, thereby embracing new situations that 
have not been considered up to now. 
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FIG. 1. (a) SNR corresponding to f(x,t) given through 
Eq. (^) with parameter values ki = 100, k% = 1, T = 4, 
j3 = 1, a = 0.5 and u>o/2n — 1. (b) Signal (filled circles) and 
noise (empty squares) for the previous situation, (c) SNR as 
in Fig. (a) but with T = 10. 
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FIG. 2. SNR corresponding to f(x,t) given through Eq. 
(^) parameter values ki = 100, = 1, with parameter values 
k = 100, T = 4, [3 = 201, a = 100 and lu /2tt = 1. 
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FIG. 3. (a) Potential Vo(x) corresponding to the force 
— fo(x)x with parameter values fci = 20, &2 = —80, — 20, 
fc 4 = -80, fc 5 = 0, ci = 1, c 2 = 1.73, c 3 = 3.16, and c 4 = 4.47. 
(b) SNR corresponding to f(x,i) given through Eq. (|Io|). 
The parameter values for fo(x) are the same as in case (a). 
Moreover (3 = 121, a — 100, and u)q/2k — 1. (c) Same 
situation as in case (b), except fcs = 20. (d) Same situation 
as in case (b), except c 4 = 7.75. 
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